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Resumen 


El presente artículo describe un enfoque original para determinar curvas 
de fragilidad de bordos perimetrales, ante la falla por erosión debida al 
desbordamiento basado en la simulación de variables aleatorias. El 
análisis de fragilidad toma en cuenta la incertidumbre asociada con la 
estimación de los parámetros de un modelo que describe la erosión 
superficial del material del bordo, y después, las curvas de fragilidad se 
expresan en función del tirante de agua, por encima de la cresta del 
bordo. En los análisis de fragilidad, tradicionalmente se supone a priori 
una distribución lognormal para las curvas de fragilidad y sus 
parámetros se obtienen posteriormente, utilizando alguna técnica de 
estimación estadística. Los resultados del presente artículo muestran 
que, al considerar una densidad de probabilidad no gaussiana para el 
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coeficiente de erodibilidad, las curvas de fragilidad difieren de manera 
significativa de la distribución lognormal. 


Palabras clave: curvas de fragilidad, erosión superficial, simulación de 
variables aleatorias y coeficiente de erodibilidad. 


Abstract 


This paper describes an original approach to determine fragility curves 
for the surge only overtopping erosion failure of levees on the basis of 
the simulation of random variables. The fragility analysis accounts for 
uncertainty in estimating the parameters that describe surficial erosion 
of the levee material and then, fragility curves are expressed as a 
function of the height of water above the levee crest. Common practice 
assumes a prior lognormal distribution for the fragility curves and their 
parameters are then obtained by some statistical estimation technique. 
The results of this paper shown that when a non-Gaussian probability 
density function is considered for the erodibility coefficient, the fragility 
curves differ significantly from the lognormal distribution. 


Keywords: fragility curves, surficial erosion, simulation of random 
variables, erodibility coefficient. 
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Introducción 


La erosión que causa el agua al fluir sobre la superficie de un bordo 
perimetral puede provocar su rompimiento, inundando súbitamente a la 
población, agricultura, ganadería e infraestructura que se desea 
proteger. En la literatura se han documentado varios casos de 
rompimiento de bordos por desbordamiento (Conagua, 2011; Cenapred, 
2014). 


Las curvas de fragilidad se han utilizado extensamente en la evaluación 
probabilista del riesgo estructural desde finales de la década de los 
setentas (Kennedy, Cornell, Campbell, Kaplan, 8€ Perla, 1980). 
Posteriormente se utilizaron en el diseño de elementos de cimentación 
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(Popescu, Deodatis, £ Nobahar, 2005). Actualmente es común analizar 
la fragilidad de bordos para protección de inundaciones (Vorogushyng, 
Merz € Apel, 2009; Schweckendiek € Vrouwenvelder, 2014 y 
Bachmann, Huber, Johann, €  Schúttrumpf, 2013). Diferentes 
organismos internacionales incluyen en sus recomendaciones el uso de 
este concepto en la evaluación probabilista del riesgo (USACE, 2010; 
SPANCOLD, 2013; CIRIA, 2013). 


Definición de curva de fragilidad 


Sea x'=[X,,X,,K,X,] un vector de estado que contiene variables 
aleatorias básicas de diseño X,,X,,K,X, con densidades de probabilidad 
conjunta N-variadas, representadas por f(x). Por simplicidad, los 


eventos de seguridad y falla ante un mecanismo específico se describen 
en términos de una función de estado g(x), escrita de tal manera que: 


leGo<¿), representa el evento de seguridad; por ejemplo, (xeS) y 
leGo >¿), representa el evento de falla; por ejemplo, (xe 3); donde ¿£ 
es un umbral que delimita las regiones de seguridad y falla. 


La probabilidad de falla dado un valor específico h de la intensidad es la 
fragilidad y se puede expresar como (Casciati 8 Faravelli, 1991): 


F,,(n)= Plx e SH = h]= 7,6) (1) 
3|n 
La curva dada por F,(h) es la curva de fragilidad y expresa la 


probabilidad de que la capacidad sea menor que la demanda dada una 
intensidad h de la solicitación. El dominio de falla depende de la 
intensidad. Conforme la intensidad aumenta, la probabilidad de falla se 
aproxima a la unidad. 


Enfoque propuesto basado en la simulación de 
variables aleatorias 
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En este contexto se define una curva de fragilidad empírica para 
diferentes valores de » en la siguiente forma: 


F, (n)= Plx e 3H =n]= F76)£, 69) (2) 


e|h 


donde 0, representa el dominio de las variables aleatorias en x e /(x) 
es una variable aleatoria binaria tal que: /(x)=0, Si g(x)<£ O I(x)=1, Si 
2(0)>6. El símbolo “4” encima de F,(h) enfatiza que la curva dada por 
la ecuación (2) es una curva de fragilidad empírica (punto por punto). 
La versión teórica de la curva empírica se puede obtener por medio de 
algún procedimiento de ajuste (Popescu et a/., 2005). 


Para determinar el valor que toma la variable aleatoria binaria es 
necesario simular las variables aleatorias del vector x. Para realizar esta 
tarea, se propone el siguiente algoritmo, basado en el método de 
descomposición de Cholesky (Rao, 1992): 


1. Especificar la matriz de correlación deseada en términos de p.. 

2. Obtener la matriz de correlación R en términos de »r, tal que: 
r =2sen(rp, /6) . 

3. Realizar la descomposición de Cholesky de kr, tal que: R=IT". 

4. Simular Ny variables aleatorias gaussianas independientes: 
Z AZ Zoro Zy) » 

5. Hacer: y=Iz. 

6. Obtener: U=0(y); donde D() es la función de distribución 
gaussiana. 

7. Hacer: X,=F (U), X,=FE,U)) y.., Xy =Fy (Uy); donde F;'() para 
j=1,2,...N son los inversos de las funciones de distribución 
deseadas. 


Con el procedimiento anterior es posible determinar la probabilidad de 
falla para diferentes valores de la intensidad h por medio de la ecuación 
(2). Note que en el paso 1 la correlación entre variables aleatorias se 
expresa en términos del coeficiente de correlación de Spearman (p,) y 


que en el paso 2 la matriz R se obtiene en términos del coeficiente de 
correlación de Pearson (r). El coeficiente de correlación de Pearson se 
utiliza para expresar la correlación entre las densidades de probabilidad 
conjuntas, que en el presente trabajo se suponen gaussianas, a partir 
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del coeficiente de correlación de Spearman, que expresa la correlación 
deseada entre las variables marginales no gaussianas (Genest € Favre, 
2007). 


Análisis de la falla por erosión debida al 
desbordamiento 


Con base en la ecuación del exceso de esfuerzo cortante (Partheniades, 
1965), la función de estado para la falla por erosión debida al 
desbordamiento de un bordo perimetral se puede escribir en la siguiente 
manera: 


23h) =E—k¿(t, -7,):1 (3) 


donde £: erosión [mm], +: tiempo [h], k,: coeficiente de erodibilidad 
[mm/h/Pa], 7.: esfuerzo cortante crítico [Pa], 7,: esfuerzo cortante en 


Cc 


la interfaz suelo-agua [Pa] (Nadal €: Hughes, 2009): 
To, = pgeh,sen0O (4) 
donde p: densidad del agua [kg/m*], g: aceleración de la gravedad 


[m/s?], 0: ángulo de inclinación del talud aguas abajo [*] y h,: tirante 
0 


normal al talud [m] (Powledge, Ralston, Miller, Chen, Clopper, 8 
Temple, 1989): 


donde »h : altura del desbordamiento [m] y n: coeficiente de Manning 
[s/m*J, 


En la ecuación (3), k, y 7, son propiedades que dependen del material 
que constituye el bordo y se determinan en pruebas de erosión de 
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campo o laboratorio; k, determina la velocidad con la que ocurre la 
erosión del bordo y 7. establece el umbral a partir del cual comienza la 
erosión del bordo. Como se ilustra en la Figura 1, cuando 7,>z., 


comienza la erosión del bordo y ésta ocurrirá a una velocidad dada por 
k,. Los valores de k, y 7, utilizados en el presente análisis de fragilidad 
se estimaron con base en los resultados de las pruebas de erosión de 
laboratorio reportadas por Briaud, Chen, Govindasamy y Storesund 
(2008) y Yuan, Tang, Li, Pan, € Amini (2015) para suelos finos 
compactados. 


ty 1)<T, 


To > To 


Figura 1. Ilustración del proceso de erosión del bordo. 


Aplicación al caso de un bordo hipotético 


En el siguiente ejemplo ilustrativo se analiza la fragilidad de un bordo 
perimetral de material fino compactado de 5 m de altura y 4 m de 
ancho, sometido a un desbordamiento en el intervalo: 0O<h<2 m. Se 


considera una situación de la práctica en la que el diseñador tiene 
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incertidumbre en la estimación de cuatro variables básicas de diseño: 0, 
n, T. Y k,. Las fragilidades se determinan para diferentes muestras de 


H,, donde H, es una variable aleatoria que refleja las fluctuaciones del 
nivel de agua en el embalse. 


Los resultados de ensayos de laboratorio en muestras de materiales 
finos compactados reportados por Briaud et al. (2008) y Yuan et al. 
(2015) se consideraron en la estimación de los valores esperados y 
dispersiones de 7, (esfuerzo cortante crítico) y k, (coeficiente de 


erodibilidad). Al construir el histograma de los valores reportados para 
7, se encontró que la densidad lognormal describe satisfactoriamente los 


resultados de los ensayes. Se observó además que el parámetro k, 


fluctúa en varios órdenes de magnitud y que su valor mínimo siempre es 
mayor que cero. Por lo tanto, se adoptó la densidad de Weibull para 
representar la incertidumbre en este parámetro. En la práctica, por lo 
general se selecciona un valor esperado para el coeficiente n de 
Manning, en función del tipo de material del bordo y se acepta que 
puede variar dentro de un intervalo conocido (Nadal € Hughes, 2009). 
La incertidumbre en este caso se modeló a través de la densidad beta 
descrita en Oboni y Bourdeau (1985). Por simplicidad, las posibles 
variaciones del ángulo e del talud aguas abajo se representaron por 
medio de una variable aleatoria gaussiana. Los parámetros del conjunto 
de variables aleatorias básicas de diseño se indican en la Tabla 1. 


Tabla 1. Parámetros probabilistas de las variables básicas de diseño. 


Parámetro Unidad EDO Cv Mín Máx. Tipo 
0 rad 0.197396| 0.1 =..- ---- |Normal 
n s/mi/? 0.0158 0.1 | 0.01 | 0.035 | Beta 
Te Pa 16 6 0.0 ---- | Lognormal 
Ka mm/h/Pa 2.5 6 0.1 ---- | Weibull 


Se ha considerado que existe correlación entre los tres parámetros que 
dependen del tipo de material del bordo (n, 7, y k,). Esta decisión no es 
subjetiva, se basa en los resultados de las pruebas de erosión de 
laboratorio reportadas por Briaud et a/. (2008) y Yuan et al. (2015). Se 
consideró una correlación negativa (-0.7) entre 7, y k, para representar 
que cuando una variable tiende a crecer, la otra tiende a disminuir, 
como se observa en los resultados de laboratorio reportados por Briaud 
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et al. (2008). Se consideró además, una correlación negativa (-0.8) 
entre 7. y n para representar que conforme el coeficiente n de Manning 


aumenta, el esfuerzo cortante crítico 7. disminuye. Esta tendencia se 


fundamenta con base en que los suelos compactados con mayor 
contenido de arena suelen exhibir coeficientes n de Manning mayores, 
pero son menos resistentes a la erosión. Los valores numéricos 
asignados a los coeficientes de correlación mencionados arriba (-0.7 y - 
0.8) reflejan el grado de dependencia que el diseñador asigna, de 
manera parcialmente subjetiva, a dichas tendencias. Como ya se 
mencionó, esta decisión se basa en el comportamiento observado de 
dichos parámetros en las pruebas de laboratorio reportadas en la 
literatura. Para ser congruentes con nuestro razonamiento, deben existir 
entonces una correlación entre n y k, pero de signo opuesto (0.8). 


Con base en las experiencias reportadas por Hughes (2010), en el 
presente trabajo se adopta una escala de erosión permisible arbitraria 
en la siguiente manera: 


e Daño menor: Si ¿>0.30 m. 
e Daño intermedio: Si £>0.60 m. 
e Daño severo: Si £>0.90 m. 


La fragilidad del bordo ante el desbordamiento se analiza para cada 
nivel de daño, especificado arriba y para dos tiempos de permanencia 
del desbordamiento; es decir, para un periodo de 4 h y otro de 20 h. El 
primer periodo representa un tiempo típico de análisis y el segundo 
periodo una condición extrema. Cada punto de la curva de fragilidad 
empírica se determinó con base en 500 000 simulaciones por medio de 
la ecuación (2). 


Resultados y discusión 


Las Figuras 2a y 2c muestran con símbolos las curvas de fragilidad 
empíricas correspondientes al primer análisis de fragilidad. Note que 
conforme la duración del desbordamiento aumenta, la probabilidad de 
falla aumenta. Se puede observar además que conforme el indicador de 
falla £ aumenta, la probabilidad de falla disminuye porque la frecuencia 
con la que ocurre una falla más severa es menor si se mantienen 
constantes los tiempos de análisis. Para un desbordamiento próximo a 2 
m de altura que permanece durante 20 h, el bordo alcanza la falla con 
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probabilidad muy cercana a uno. Siguiendo a Shinozuka (2000), las 
curvas de fragilidad en las Figuras 2a, 2b y 2c se pueden interpretar 
como la probabilidad de que el bordo sufra al menos el daño 
especificado cuando se somete a un desbordamiento con intensidad h,. 
Las probabilidades de falla aun para el periodo típico de análisis de 4 h 
son, en efecto, elevadas. Considerando que en los análisis se han 
utilizado valores representativos de las propiedades de erosión de los 
suelos compactados, por lo tanto estas probabilidades también son 
representativas de un bordo sometido al desbordamiento y pueden 
explicar las observaciones publicadas por algunos autores. 
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Figura 2. Curvas de fragilidad para diferentes condiciones de falla: a) 
Curva de fragilidad para un nivel de daño menor; b) Curva de fragilidad 
para un nivel de daño intermedio; c) Curva de fragilidad para un nivel 
de daño severo. 
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Por ejemplo, Gilbert y Miller (1985) han reportado que la erosión de 
bordos de arcilla compactada normalmente inicia entre 10 h y 20 h 
antes de su ruptura. Los resultados del presente análisis de fragilidad 
también pueden explicar la frecuencia con que ocurren fallas en presas 
por desbordamiento (48.4%) y tubificación (46.1%) (Foster, Fell, 8 
Spannagle, 2000). Recordemos que la erosión del talud aguas abajo 
favorece la tubificación de los terraplenes de materiales finos 
compactados sometidos a un flujo de agua en su interior. Para aumentar 
la confiabilidad de los bordos ante el desbordamiento, algunos autores 
muestran las ventajas de sembrar césped sobre el bordo (Powledge et 
al., 1989); o mejor aún, incrementar la resistencia a la erosión del 
césped sembrado utilizando geomalla (Yuan et al., 2015). 


Las curvas de fragilidad teóricas asociadas a las curvas empíricas se 
determinaron utilizando la técnica de mínimos cuadrados (Popescu et 
al., 2005). Se encontró que cada curva teórica está dada por una 
distribución de Weibull compuesta (Razali 8 Al-Wakeel, 2013) en la 
siguiente forma: 


Gy (0) = 2 w,F, 01) (6) 


donde G(x): función de distribución compuesta, w,: ponderación 


asociada a la ¡-ésima distribución de Weibull y F,(x): función de 
distribución de Weibull dada por: F,(x)=1-exp[-(x/a)] con parámetros 
de escala y forma dados por a y £, respectivamente. 


Las curvas de fragilidad teóricas se muestran con líneas discontinuas en 
las Figuras 2a, 2b y 2c. Los parámetros que definen estas curvas 
teóricas se indican en la Tabla 2. Note que las curvas teóricas describen 
exactamente la curva empírica. Dada la forma de la curva teórica se 
puede concluir que la densidad de probabilidad del parámetro k, 


(coeficiente de erodibilidad) ejerce la mayor influencia sobre la forma de 
la curva de fragilidad y que la influencia de n (coeficiente de Manning) y 
r. (esfuerzo cortante crítico) es menor que la influencia de k,. La forma 


de la curva teórica dada por la ecuación (6) se puede explicar por la no- 
linealidad del modelo conjunto de erosión dado por las ecuación (3), 
ecuación (4) y ecuación (5). Los cálculos realizados durante el presente 
trabajo mostraron que la correlación entre las variables de diseño 
modifica los parámetros de la curva teórica, pero no influye en su tipo. 
Algunos autores han sugerido el uso de densidades exponenciales para 
los parámetros de erosión de los suelos (Govindaraju, 1998). Otra 
posibilidad es considerar densidades de probabilidad lognormales para 
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representar la incertidumbre en la estimación de estos parámetros. Sin 
embargo, en el caso del coeficiente de erodibilidad, una densidad 
lognormal no es representativa del comportamiento físico de dicho 
parámetro, porque la erodibilidad de los materiales finos compactados 
no puede ser nula; parte de un valor mínimo mayor que cero y estará 
más alejado del cero conforme las incertidumbres en la calidad de los 
materiales del banco de préstamo y los procedimientos de colocación de 
los materiales sean mayores. Dado el efecto del tipo de densidad de k, 


sobre la forma de la curva de fragilidad, es recomendable evaluar sus 
consecuencias en un análisis probabilista del riesgo en cada caso 
particular (Tabla 2). 


Tabla 2. Parámetros descriptivos de las curvas de fragilidad. 


Ñ Dura- Ponderación Escala Forma 
Daño des 

ción Wi W2 W3 0 0 03 Pi Ba B3 

Menor 4 h 0.179 | 1.303 | 0.277 | 0.065 | 3.867 | 0.672 | 0,804 | 3.120 | 1.167 
20h | 0.184 | 0.551 | 0.263 | 0.419 | 0.112 | 0.434 | 9.029 | 0.601 | 3.014 

atemnadió 4 h 0.196 | 0.120 | 0.181 | 0.163 | 1.459 | 1.898 | 0.762 | 0.981 | 1.405 
20h | 0.414 | 0.505 | 0.208 0.893 0.255 9.815 | 5.872 | 0.636 | 0.355 

Sevsio 4h 0.083 | 0.354 | 0.037 | 0.936 | 1.776 | 0.102 | 0.613 | 0.714 | 0.714 
20h |0.345 | 0.525 | 0.446 1.392 3.560 0.686 | 7.480 | 0.571 | 0.516 

Conclusiones 


El presente artículo mostró que las curvas de fragilidad para la falla por 
erosión de un bordo perimetral sometido a un desbordamiento difieren 
significativamente del modelo lognormal cuando se utiliza un modelo 
probabilista no gaussiano, para el coeficiente de erodibilidad del material 
del bordo. Este resultado tendrá implicaciones prácticas importantes 
cuando la curva de fragilidad se combine con la curva asociada con la 
amenaza específica y las consecuencias para cuantificar el riesgo. Por lo 
tanto, es recomendable que el efecto de las incertidumbres no 
gaussianas sobre el riesgo se evalúe en todos los casos. Este trabajo 
presentó herramientas para modelar apropiadamente este tipo de 
incertidumbres y se ponen a disposición de los interesados contactando 
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al primer autor. Finalmente, se recomienda a los organismos nacionales 
con responsabilidad en proyectos de obras de protección contra 
inundaciones utilizar el concepto de curvas de fragilidad en la evaluación 
del riesgo de sus obras de protección y encauzamiento, como sucede en 
la actualidad con diversos organismos internacionales (USACE, 2010; 
SPANCOLD, 2013; CIRIA, 2013). 
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